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Abstract. H2 column densities towards CO clouds in the 
LMC and SMC are estimated from their far-infrared sur- 
face brightness and HI column density. The newly derived 
H2 column densities imply N(H2) / I(CO) conversion fac- 
tors (in units of 10^^ niol cm^'^ (Kkms^^)^^) Almc = 
1.3±0.2 and Asmc = 12±2. LMC and SMC contain total 
(warm) H2 masses of l.OiO.S x 10^ Mq and 0.75±0.25 
X 10* Mq respectively. Local H2/HI mass ratios simi- 
lar to those in LMC and SMC are found in the magel- 
lanic irregulars NGC 55, 1569, 4214, 4449 and 6822 and 
in the extragalactic HII region complexes NGC 604, 595 
and 5461 in M 33 and M 101 respectively. In these HII 
regions and in NGC 4449, we find X= 1-2; in NGC 55, 
4214 and 6822 X = 3-6 again in units of lO^^ mol cm"^ 
(Kkms^^)^^. The post-starburst galaxy NGC 1569 has a 
very high value similar to that of the SMC. 

The CO-H2 conversion factor X is found to depend 
on both the ambient radiation field intensity per nucleon 
ctfir/ and metallicity [0]/[H]: log Aoc 0.9±0.1 log-^ 
- 3.5±0.2 log|^. Neglecting dependency on radiation field, 
a reasonable approximation is also provided by log A oc - 
2.7±0.3 logjgj. Milky Way values are consistent with these 
relations. This result is interpreted as the consequence of 
selective photodissociation of CO subjected to high ra- 
diation field energy densities and poor (self)shielding in 
low-metallicity environments, and especially the preferen- 
tial destruction of diffuse CO in 'interclump' gas. 

Although locally H2 may be the dominant ISM- 
component, the average global H2/HI mass ratio is 
0.2±0.04 and the average H2 gas mass fraction is 
0.12±0.02. Magellanic irregulars have warm molecular gas 
fractions very similar to those of our Galaxy, whereas 
other global properties (mass, luminosity, metallicity, CO 
luminosity) are very different. 

Key words: Galaxies: individual: LMC; SMC - Galaxies: 
ISM; irregular; Magellanic Clouds - Infrared: ISM: contin- 
uum - ISM: molecules 



1. Introduction 

1.1. H2 content of galaxies 

The existence of molecular hydrogen (H2) in interstellar 
space was suggested as early as 60 years ago by Eddington 
(1937) and Stromgren (1939). Thirty years later, Gould 
& Salpeter (1963) and HoUenbach et al. (1971) predicted 
that it could be a large fraction of all hydrogen. However, 
H2 is difficult to observe directly, because it is a symmet- 
rical molecule lacking a dipole moment. Nevertheless, it 
has been observed in absorption at UV wavelengths and 
in emission at infrared wavelengths. Because the emission 
arises mostly in warm or hot molecular gas, it has been 
virtually impossible to deduce total amounts of H2 which 
is expected to be present mostly at low temperatures. 

As H2 is an abundant and important constituent of the 
interstellar medium in galaxies, there has been great inter- 
est in determining its presence and properties. Because it 
is commonly assumed that star formation requires inter- 
stellar clouds to pass through a cool, high-density phase in 
which most of the hydrogen is in molecular form, studies of 
star formation in external galaxies also seek to determine 
H2 amounts in such galaxies. 

Emission from the tracer CO molecule has been and 
still is widely used to determine the distribution and 
amount of molecular hydrogen in our own Milky Way and 
other galaxies. Although the usually optically thick -'^^CO 
emission does not provide direct information on column 
densities, empirical relations between CO luminosities and 
virial masses of molecular clouds in the Galaxy suggest 
that circumstances nevertheless allow its use in an indirect 
manner. The underlying thought is that CO is distributed 
in a very clumpy manner, and that clumps are not self- 
shadowing. The strength of the signal received from many 
clumps in a single observing beam thus provides a mea- 
sure for their total projected area weighted by brightness 
temperature. If the clumps are statistically similar from 
one line of sight to the other, we thus have a measure for 
the number of clumps per beam area, hence for the to- 
tal amount of molecular material. In fact, the high optical 
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depth of ^^CO emission is then a boon, as it makes such 
determinations to first approximation independent of the 
actual [C0]/[H2] abundance. 

1.2. Problems with CO-based methods 

The most commonly used methods to estimate molecu- 
lar hydrogen contents of extragalactic objects are cither 
application of a 'standard' CO to H2 conversion factor X 
(defined as the ratio of molecular hydrogen column density 
N(H2) to velocity-integrated CO intensity I(COJ) derived 
from Milky Way observations, or application of the virial 
theorem to observed CO clouds. The first method assumes 
similarity of extragalactic molecular clouds and Galactic 
clouds, or at least that the effects of different physical con- 
ditions cancel one another. In environments that are very 
different from those in the Galaxy, such as those found in 
galaxy central regions or in low-metallicity dwarf galaxies, 
this method must be considered unreliable (c/. Elmegreen 
et al. 1980; Israel 1988; Maloney & Black 1988; Elmegreen 
1989; Maloney 1990a). For instance, application of this 
method to the very low CO luminosities commonly ob- 
served for irregular dwarf galaxies would suggest negligi- 
ble amounts of H2 (Israel et al. 1995) and consequently 
unusually high star formation efficiencies (Israel 1997). In 
contrast, direct evidence for X factors varying by more 
than an order of magnitude, probably as a function of lo- 
cal conditions, has been presented for Galactic clouds by 
Magnani & Onello (1995). We thus agree with Roberts & 
Haynes (1994): 'values of the molecular hydrogen content 
in late-type systems derived in this manner are uncertain 
and possibly too low by up to an order of magnitude' 

The second method frequently used estimates total 
molecular cloud mass from observed parameters such as 
CO extent R(CO) and velocity dispersion <:\v(CO). Al- 
though this method, not assuming similarity between 
Galactic and extragalactic clouds, is preferable, it is like- 
wise beset by problems, as it requires correct determina- 
tion of the structure and dynamics of the observed clouds. 
For instance, the value of the virial constant used to con- 
vert observed parameters into mass may vary by a fac- 
tor of four depending on the assumed condition of the 
system (see e.g. McLaren et al. 1988; McKee & Zweibel 
1992), while it is not clear that the virial theorem is in 
fact relevant. If one considers the morphology of molecular 
complexes such as the ones in Orion (Bally et al. 1987), 
Taurus (Ungcrechts & Thaddeus 1987) or indeed in the 
LMC (Israel & de Graauw 1991; Kutner et al. 1997) it 
is hard to imagine that these very elongated structures 
with little systematic velocity structure actually represent 
virialized complexes. Maloney (1990b) has shown that the 
correlation between CO luminosities and virial masses of 
Galactic molecular clouds follows directly from the size- 
linewidth relationship exhibited by molecular clouds and 
does not require virial equilibrium at all. Molecular hydro- 
gen masses have also been determined applying X-factors 



scaled from Xcai by L( CO) as a function of Av (e.g. LMC 

- Cohen et al. 1988; SMC - Rubio et al. 1991). 

Especially in the large linear beamsizes typical of 
extragalactic observations, actually unrelated clouds at 
somewhat different velocities may blend together, lead- 
ing to unrealistical values of both cloud complex radius 
R and velocity dispersion dv. The derived (virial) masses 
may then either overestimate or underestimate the actual 
mass, depending on circumstances. For instance, consider 
an area mapped with a large beam blending together N 
unrelated clouds, each having a true mass Mvir = a r dv^. 
Here, Tq is the diameter of a single cloud and dvo its ve- 
locity dispersion. The true total mass is thus N aro dv^. 
Cloud emission is measured over an area with radius R 
= iV^/^ b To in which b Tq is the projected separation be- 
tween cloud centers. Unjustified application of the virial 
theorem on this observation suggests a total mass Mvir = 
a iV^/^ b To dvf, where dvi now is the dispersion derived 
from the velocity width of the sum profile of all clouds 
within radius R. The ratio of the derived mass over the 
true mass is thus: 

Mderivcd/Mtrue = b N-^/^ (dvi/dVo)^ (1) 

U N < b'^ and dt;i > dVo, this will result in a poten- 
tially large overestimate of the mass. However, if instead 
the unrelated clouds are at more or less identical radial ve- 
locities, dvi « d^^o, the true mass is underestimated if N 
> b^. Such a situation may occur in low-metallicity dwarf 
galaxies with relatively small velocity gradients. It occurs 
if we have a large number of clouds with small projected 
distances; a more physical equivalent is a very filamentary 
structure of the molecular material. 

A further problem in estimating H2 masses from CO 
observations is the need to assume virtually identical dis- 
tributions for both. If CO is significantly depleted, H2 may 
well occur outside the area delineated by CO emission and 
its amount is underestimated by the CO measurements. 
This effect appears to lie at the base of the size dependence 
of N(H2)/I(C0) ratio, noted by Rubio et al. (1993) and 
Verter & Hodge (1995). In low-metallicity galaxies suffer- 
ing CO depletion, this results in a lack of CO emission 
in complexes observed on large angular scales. Observa- 
tions on small angular scales selectively concentrate on 
the densest molecular components, that have resisted CO 
depletion most effectively, so that the N(H2)/I(C0) ra- 
tio looks progressively more 'normal' notwithstanding the 
lack of CO in most of the complex. 

Thus, in order to estimate H2 content of such galax- 
ies, it is desirable to use a method that does not require 
specific assumptions on or knowledge of the; detailed struc- 
ture and dynamics of the molecular clouds involved. Use 
of far-infrared data in principle provides such a method 
(e.g. Thronson et al. 1987, 1988; Israel 1997). 

2. Method and data 
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2.1. Estimating N(H2) from ctfir and N(HI) 

H2 column densities are derived in the manner used on 
NGC 6822 by Israel 1997. At locations well away from 
star-forming regions and CO clouds, the ratio of neutral 
hydrogen column density to far-infrared surface brightness 

{N(HI)/a-piii)o is determined. In the absence of molecu- 
lar gas, {N(HI)/aFi^)o equals iVn/cTFiR, which is a mea- 
sure for the ambient gas-to-dust ratio. The observed ctfir 
values at locations that contain H2, as betrayed by CO 
emission, reduced to the reference dust temperature Tdo 
and then multiplied by (N(HI)/ (7-Fm)o thus provide the 
total hydrogen column density A^h • The actual gas-to-dust 
ratio, which depends on poorly known dust particle prop- 
erties, does not need to be known as long as it does not 
change from source to reference position. In the small ir- 
regular galaxies considered here, abundance gradients are 
negligible (cf. Vila-Costas & Edmunds 1992), so that we 
may safely assume no change in gas-to-dust ratio as a func- 
tion of position in the galaxy. When A'h is known, N(H2) 
is found by subtracting the local N(HI) value: 

2 N(H2) = [{N(HI)/ aFiR)o m ctfir] - N(HI) (2) 

In eqn. 2, f(T) is a function which corrects ctfir for 

the emissivity difference due to the (generally small) dif- 
ference of Td from Tdo; for small temperature differences 
f(T) is close to ( Tdo/ T^)^. Temperatures Td are derived 
from the IRAS 60/im/lOO/im flux ratio assuming a wave- 
length dependence for emission Ix ex A~"iJA- Here and 
in the following we will assume n = 2. The temperature 
correction assumes that the number distribution of dust 
particles emitting at varying temperatures does not differ 
significantly from one location to another. This is a reason- 
able assumption for values of J( T) not too far from unity, 
but may introduce significant errors for very large or very 
small values oif(T). The CO to II2 conversion factor X fol- 
lows from the observed CO strength: X = N(H2)/I(C0). 

This method of estimating H2 column densities de- 
pends on observed quantities independent of the actual 
spatial or kinematical distribution of the molecular m,ate- 
rial. It has this property in common with the methods 
used by Bloemen et al. (1986) and Bloemen et al. (1990) 
to estimate the same quantities in the Milky Way galaxy. 
It avoids the major weakness of the virial method dis- 
cussed above, as there is no need to determine the struc- 
ture of the molecular cloud complexes, to separate unre- 
lated clouds in the same line of sight, or even to resolve 
the molecular clouds. It is important to emphasize that in 
this method, the absolute gas-to-dust ratio plays no role, 
nor does the the actual dust mass. We thus avoid a major 
uncertainty associated with other infrared-derived H2 es- 
timates, where the infrared flux is used to calculate a dust 
mass, which is then converted into a gas mass. Likewise, 
our results are independent of CO measurements, and as 
we will show below, the observational uncertainties are no 



worse than those associated with the traditional methods 
and probably better. 

The column densities N(H2), and consequently X, de- 
termined in this paper are properly lower limits (Israel 
1997). (i). If some II2 were to be present at the null posi- 
tions where we assumed none, the total hydrogen column 
density corresponding to unit infrared luminosity is un- 
derestimated, implying higher actual N(H2) values than 
derived, (ii). If, unexpectedly, the /lotter infrared sources 
were to be relatively rich in cooler dust, the observed in- 
frared emission does not sample the total amount of gas, 
hence N(H2), will be higher than estimated. (Hi). If, in 
regions of bright infrared emission, higher radiation den- 
sities would cause increased dust depletion, these regions 
will be characterized by a higher gas-to-dust ratio than as- 
sumed, again leading to higher than derived actual N(H2) 
values. This is expected only to be important for HII re- 
gions filling a significant fraction of the beam. 

Errors in the assumptions would thus cause N(H2) and 
Xto be higher rather than lower. Although these errors are 
hard to quantify, we consider it unlikely that their effect 
will exceed a factor of two. The calculated total hydro- 
gen column densities iVfj carry with them the combined 
uncertainty in the determinations of {N(HI)/ apinjo, f(T) 
and (Tfir- Because these quantities are compared in a rel- 
ative rather than an absolute sense, the uncertainty A iVu 
is of the order of 20% - 30% for the cases discussed be- 
low. The uncertainty in the calculated values of N(H2) is 
larger. Since the N(HI) determinations are considered to 
be rather accurate, it depends on the molecular to atomic 
hydrogen ratio: A7V(i?2) = A ArH(l + 0.5 N(HI)/N(H2)) 

Thus, for II2 column densities equal to or higher than 
those observed in HI, the relative H2 uncertainty is typi- 
cally less than 50%. For HI column densities substantially 
higher than the derived H2 column density, the relative 
uncertainty may become considerable. However, this situ- 
ation almost exclusively occiirs at low absolute N( H2) val- 
ues where a relatively large uncertainty still corresponds 
to an acceptable uncertainty in the absolute value. The 
uncertainty in the derived value of X, in turn, includes 
both the uncertainty in N(H2) and in I (CO). Since the 
latter is usually much smaller than the former, the uncer- 
tainty in X is actually dominated by that in N(H2). The 
combined effect of uncertainties in the observational values 
and in the assumptions implies a rough overall uncertainty 
of about a factor of two for individual determinations. 

2.2. Data and results 

All data were taken from the literature or existing 
databases. The CO, HI and far-infrared data included in 
the comparison are selected to have similar resolutions. 
This resolution is determined by the lowermost resolution 
to which the other data are degraded, if necessary. 
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Table 1. LMC Data (unit area 0.043 kpc^) 
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6.0 


0.38 


0.7 


2.6±0.7 


7.4 


0.65 


0.75 


3.4±1.1 


10.3 


N83B 


6 


12 


0.26 


1.5 


2.4±0.8 


3.2 


0.56 


1.15 


2.1±0.9 


17.2 


Nil 


6 


1.2 


1.00 


1.2 


0.8±0.3 


1.3 


0.42 


1.35 


0.6±0.2 


4.4 


Nil-North 


11 


1.2 


1.00 


1.1 


0.8±0.3 


1.5 


0.44 


0.75 


1.1±0.7 


4.5 


(Bar) 


13 


4.8 


0.38 


1.3 


1.4±0.5 


2.1 


0.51 


0.75 


1.9±0.9 


11.7 


N105A 


15 


2.4 


0.57 


1.0 


1.0±0.4 


2.0 


0.49 


1.15 


1.0±0.5 


7.8 


N113 


- 


1.8 


0.57 


1.2 


0.6±0.3 


1.0 


0.37 


0.40 


1.4±1.1 


7.8 


N30 


18 


2.1 


1.00 


0.9 


1.9±0.6 


4.2 


0.60 


1.15 


1.6±0.6 


4.5 


(Bar) 


19 


12 


0.31 


2.0 


3.2±1.0 


3.2 


0.56 


2.50 


1.3±0.4 


14.3 


N44 


20 


1.8 


0.51 


1.5 


0.3±0.3 


0.4 


0.21 


1.15 


0.3±0.3 


8.7 


N18/N144 


22 


1.2 


0.57 


0.7 


0.4±0.2 


1.1 


0.40 


0.75 


0.6±0.3 


7.8 


N132 


23 


3.6 


0.38 


1.5 


0.8±0.4 


1.1 


0.38 


1.50 


0.5±0.3 


11.7 


N48 


26 


4.2 


0.38 


1.0 


1.3±0.4 


2.6 


0.53 


0.40 


3.3±1.8 


11.7 


N206 


27 


2.4 


0.51 


1.6 


0.6±0.3 


0.8 


0.32 


1.15 


0.5±0.3 


8.7 


N148/N150 


29 


4.5 


0.38 


1.1 


1.4±0.5 


2.6 


0.53 


0.40 


3.6±1.9 


11.7 


N57 




i.O 


U.O/ 


1 n 
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U.4o 


n Ac\ 
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'7 Q 
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IN 04 


32 
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0.08 


3.2 


9.7±2.7 


6.1 


0.64 


1.15 


8.4±2.9 


52.9 


30 Dor 


32'^ 


9.0 


0.26 


2.5 


1.4±0.6 


1.1 


0.39 


1.50 


0.9±0.4 


17.1 


N157B 


33^* 


12 


0.26 


3.2 


1.9±0.8 


1.2 


0.40 


4.00 


0.5±0.3 


17.1 


N159 


34 


1.2 


1.20 


3.0 


0.1±0.4 


0.1 


0.05 


0.55 


0.2±0.6 


1.9 


N214 


35 


3.6 


1.00 


3.5 


2.3±1.0 


1.3 


0.44 


4.50 


0.5db0.2 


5.1 


N176 


36 


1.8 


1.45 


3.0 


1.4±0.7 


0.9 


0.36 


2.30 


0.6±0.3 


3.1 


N216 


37 


4.8 


0.57 


3.0 


1.6±0.7 


1.1 


0.38 


2.60 


0.6±0.3 


7.8 


N167 


Mean'' 


4.3 


0.61 


1.7 


1.3±0.2 


1.9 


0.42 


1.40 


1.3±0.2 


9.2±1.0 





Notes: a. Uncertainty in I(CO) is 0.24 Kkms"^ (Cohen ot al. 1988); b. Uncertainty in is 25%; c. Henize 1956; d. Data 

unreliable because of location on strong emission gradients, e. Mean values excluding 30 Doradus. 



2.2.1. LMC 

The far-infrared data are from Schwering 1988, who con- 
veniently produced maps of infrared luminosity over HI 
mass at 15' resolution (corresponding to 235 pc) and dust 
temperature at 8' resolution (Fig. 1). The HI data (res- 
olution 15') are from Rohlfs et al. 1984. The average of 
six positions in the main body of the LMC, well away 
from CO clouds and bright HII regions is {N(HI)/ (7fir)o 
= 2.25x10^'^ cm~^/ Wm~^sr~^ (corresponding to L/ M = 
1.7 Lq/Mq) at a reference temperature T^o = 25.5 K. 
From the internal variation, we estimate its uncertainty 
to be about 10%. The uncertainty in f(T) is about 20% 
and that in itfir about 10%. 

In Table 1 we have listed data for several of the CO 
cloud complexes detected by Cohen et al. (1988) convolved 
to a resolution of 15' (e.g. Meinert, 1992). Except for cloud 
31, all CO clouds considered have diameters larger than 
15'). Weaker CO sources are included only if identifica- 
tion with an HII region complex support their validity. In 
Table 1, the first column identifies the CO cloud by its 



number in Table 1 of Cohen et al. (1988). Column 2 lists 

the far-infrared surface brightness at the peak CO posi- 
tion, column 3 the value of f(T) based on the dust tem- 
perature derived from the Fqq /Fioo flux ratio and column 
4 the HI column density. Column 5 gives the molecular 
hydrogen column densities calculated according to eqn. 2, 
while columns 6 and 7 give the resulting ratios of molecu- 
lar hydrogen to atomic hydrogen and total gas (including 
helium) respectively. Column 8 gives the integrated CO 
intensity and column 9 the resulting value of X. In col- 
umn 10 we give the ratio of the observed infrared surface 
brightness {not reduced in temperature) over total hydro- 
gen column density A^h = 2N(H2) + N(HI). This ratio is 
a measure of the ambient radiation field strength per H 
nucleon. Finally, column 11 lists HII region(s) associated 
with the molecular cloud. In most cases the HII region 
extent is much less than the 15' scale relevant to the data 
used. 

Some further comments are in order. Clouds 34, 35 
and 36 are located south of the bright HII regions as- 
sociated with the Doradus complex. Major CO emission 
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Fig. 1. LMC. Top: ratio of far-infrarod luminosity to neutral 
hydrogen column density at 15' resolution. Contours are at 
5.3, 13, 26, 52, 105 x 10"^* W m"^ sr"^ cm^. Bottom: dust 
temperature Td at 8' resolution. Contours are at 23, 28, 32 
and 36 K. 

occurs with little or no optical counterpart. There is rel- 
atively strong HI emission, but the far-infrared surface 
brightness decreases smoothly. The results for N157B and 
N159 (clouds 32 and 33) are uncertain, as both are at steep 
far-infrared gradients. N159 is also on a steep CO emis- 
sion gradient in the opposite direction. Consequently, the 
resulting value of N{H2 depends critically on the precise 
(within a fraction of the resolution) position used. The 
mean value of the CO to H2 conversion ratio (excluding 
30 Doradus) is X= 13(±2)x lO^^cm-^. Its uncertainty is 
determined by that in {N(HI) / (7Ym)o; which docs not de- 
crease with increasing sample size, whereas all other errors 
do. 

2.2.2. SMC 

The far-infrared data are from Schwering's (1988) 
maps of dust temperature and of far-infrared lumi- 
nosity over HI mass (Fig. 2). The HI data at the 
same resolution are from McGee & Newton (1981). 



Fig. 2. SMC. Top: ratio of far-infrared luminosity to neutral 
hydrogen column density at 15' resolution. Contours are at 2.6, 
5.2, 10.5, 15.7, 21, 26 x 10"^^ W m^^ sr"^ cm^ Bottom: dust 
temperature Ta at 8' resolution. Contours are at 28, 32 and 
36 K. 

For the SMC we find an average {N(HI)/ aFm)o = 
1.65(±0.25)xl0^* cm~^/ Wm"^sr^^ (corresponding to 
L/M ~ 0.23 Lq/Mq) for various positions in and near 
the bar, at a reference temperature Tdo = 28 K. In Table 
2 we list the SMC data for several of the CO cloud com- 
plexes detected by Rubio et al. (1991) in the same format 
as Table 1. 

The mean value of the CO to H2 conversion ratio is 
X = 120(±30)x 1020cm-2. The uncertainty in the null 
determination again dominates, but less decisively because 
of the relatively small sample size in Table 2. 

2.2.3. NGC 55, NGC 1569, NGC 4214 and NGC 4449 

Four other irregular galaxies have far-infrared, HI and CO 
data at similar resolutions (Table 3). For these galaxies, 
we used far- infrared data at a resolution of 1.4' obtained 
with IRAS CPC instrument at SO/xm and lOO/xm (F. Sloff, 
unpublished; Van Driel et al. 1993). For consistency, we 
interpolated the CPC 50/um fluxes to 60/im; as the abso- 
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Table 2. SMC Data (unit area 0.061 kpc^) 
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0.74 


7.8 


6.3±2.8 


1.6 


0.46 


0.42 


15±7 


0.7 


N16/N17 


SW2 


1.3 


0.67 


8.8 


3.6±2.1 


0.8 


0.33 


0.36 


10±6 


0.8 


N25/N26 




1.5 


0.55 


7.0 


4.2±2.1 


1.2 


0.40 


0.20 


20±12 


1.0 


N32/N37 


NEl 


0.3 


0.90 


3.0 


1.3±0.8 


0.9 


0.34 


0.20 


7±4 


0.6 


N72 


NE2 


1.3 


0.55 


4.3 


4.4±1.8 


2.1 


0.50 


0.33 


13±6 


1.0 


N66 


NE3 


0.9 


0.67 


5.9 


2.5±1.5 


0.9 


0.34 


0.30 


8±5 


0.8 


N76 


Mean 


1.1 


0.68 


6.1 


3.7±0.8 


1.3 


0.40 


0.30 


12±2 


0.8±0.1 





Notes: a. Uncertainty in 1(00) is 0.06 Kkms"^ (Rubio et al. 1991); b. Uncertainty in -^S is 25%; c. Henize 1956. 
Table 3. Other Galaxies 



Galaxy apia f(T) N(HI) N(H,) H™ i(CO)'^ X ^ ^ Unit 

lO"*^ 10^^ 10^^ 10"^* Area 





Wm-^sr-i 






cm ^ 






Kkms"^ 


cm ^/ Kkms ^ 


Wm ^sr ^cm^ 


kpc^ 


NGC 55 


20 


1.0 


8.0 


6.0±3.0 


1.5 


0.44 


2.0 


3.0±1.7 


10±2 


0.6 


NGC 1569 


12 


1.5 


3.5 


6.0±3.0 


3.4 


0.57 


0.37 


16±8 


8±2 


0.7 


NGC 4214 


7 


0.9 


1.6 


2.3±0.9 


2.9 


0.55 


0.7 


3.3±1.5 


12±2 


4.2 


NGC 4449 


1.6 


0.45 


1.85 


0.7±0.3 


0.6 


0.28 


0.95 


8.0±3.0 


2.3±0.6 


2.1 


NGC 6822-HII 


0.85 


0.7 


1.8 


2.8±0.9 


3.1 


0.56 


0.50 


5.5±1.1 


1.1±0.2 


0.1 


NGC 6822-IR 


0.75 


0.8 


1.5 


3.0±0.5 


4.0 


0.59 


0.85 


4.7±1.3 


1.0±0.2 


0.1 


NGC 604 


2.0 


0.6 


3.0 


2.6±1.0 


1.7 


0.47 


1.2 


2.2±0.9 


2.4±0.5 


0.5 


NGC 595 


1.1 


1.0 


1.9 


1.0±0.4 


1.1 


0.38 


0.9 


1.2±0.9 


2.7±0.6 


0.5 


NGC 5461 


6 


1.0 


1.6 


2.7±0.9 


3.4 


0.57 


2.5 


1.1±0.4 


8±2 


3.4 



Note: a. For CO details, see text. 



lute calibration of the CPC is unreliable (Van Driel et al. 
1993), we scaled all CPC fluxes by the IRAS survey fluxes. 
In the case of NGC 55 we verifled the outcome by com- 
parison with the IRAS survey image-sharpening (PME) 
result published by Bontckoe et al. (1994). 

NGC 55 was sampled at the CO cloud detected by 
Dettmar & Heithausen (1989) and at two reference po- 
sitions 1.5' on either side of this peak. Using HI data 
from Hummel et al. (1986), we found for the reference 
value {N(HI)/aFm)o = 1.4±0.2xl028 cm'V Wm-^sr"!. 
Dettmar & Heithausen (1989) give a CO surface bright- 
ness of 3 Kkms~^ and a source size of 100" x 40", so that 
I(CO) = 2 Kkms-i in a 1' beam. 

NGC 1569 was observed in CO by Greve et al. (1996) 
who detected a cloud with a peak I(CO) = 2.9 Kkms"-*^ 
and a size of 22", corresponding to a CO intensity of 
0.37 Kkms~^ in a 1' beam. HI data at 1' resolution are 
from J. Stil (private communication; see also Israel & Van 



Dricl 1990). At the reference position 1.2' to the southeast, 
where Greve et al. (1996) did not flnd CO emission, we de- 
termined {N(HI)/aFm)o = 7.7x10^8 cm'V Wm-^sr"!. 
Infrared emission gradients render this result uncertain 
by 25%. As weak CO emission might be present outside 
the limited area mapped, the uncertainty in X may be as 
high as 65%. If we take the Young et al. (1984) results 
{I(CO) = 1.1±0.3 Kkms-i in a 50" beam), we would flnd 
an X-ratio only half the value in Table 3, which we take 
as indicative of the uncertainty in X. 

NGC 4214 was observed in CO by Becker et al. (1995) 
who found a cloud complex of dimensions 36"x22" with 
a peak I(CO) = 2.8 Kkms~^ in a 21" beam. After cor- 
rection for beamsize, this is consistent with earlier and 
less accurate determinations by Tacconi & Young (1985) 
and Thronson et al. (1988). Weak CO was also detected 
by Ohta et al. (1993) in a 15" beam towards positions 
45" southeast and 35" northwest of the reference position 
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given by Becker et al. (1995). On the basis of all avail- 
able data, wc take I(CO) = 0.7(±0.15) Kkms"! in a 1' 
beam. HI data at 1' resolution are from AUsop (1979). At 
two reference positions, we determined {N(HI)/ a-Fm)o = 
1.2±0.2xl027 cm-VWm-2sr-i. 

NGG 4449 was observed in CO by Hunter & Thron- 
son 1996 (65" beam) and by Sasaki et al. 1990 (15" 
beam). Taking into account beam sizes and efficiencies, 
the data agree well. HI data at 1' resolution are from 
the WHISP database (J. Kamphuis, private communica- 
tion). We obtained reasonably accurate infrared surface 
brightnesses for Hunter & Thronson's regions 1 through 
4 only and determined {N(HI)/ (Jyik)o = 4.35±1.0xl02'^ 
cm~^/ Wm~^sr^^. Table 3 gives the mean of the individ- 
ual results for the four positions, weighted by I(CO). 

2.2.4. Extragalactic HII regions 

In Table 3, we have also included the results obtained for 
NGC 6822 (Israel 1997). Two sets of entries are given: 
NGC 6822-HII represents the mean values towards the 
HII region complexes Hubble I/III, V and X, and NGC 
6822-IR those towards the infrared sources 4 and 6 not 
associated with major HII regions. 

We also included data for the bright HII regions NGC 
604 and NGC 595 in M 33. The infrared data were taken 
from Rice et al. (1990), HI data at the same resolution 
from Deul (1988) and CO data from Blitz (1985). Note 
that the results for NGC 595 are rather uncertain, as a 
reliable flux at 100/zm is hard to determine; we assumed 
essentially similar infrared flux distribtions for both NGC 
604 and NGC 595. Finally, we also included NGC 5461, 
the brightest HII region in the galaxy M 101. Infrared 
data were taken from Bontekoe et al. (1994), HI data from 
van der Hulst & Sancisi (1988) and CO data from Blitz 
et al. (1981) and Kenney et al. (1991). In this case, we 
applied a correction factor of 1.3 to allow for the presence 
of significant amounts of HII; this effect is negligible for 
the M 33 objects. Because the parent galaxies M 33 and 
M 101 have radial abundance gradients, we selected null 
positions adjacent to the HII regions used. 

3. Analysis & Discussion 

3.1. The CO to H2 conversion factor 

In the sample galaxies, we find CO to H2 conversion factor 

X values much higher than the range 0.2-0.4 x lO^'^ cm~^ 
(Kkms~^)~^ commonly adopted for Milky Way objects. 
In the LMC, we find a mean value X = 1.3 x 10^^ cm~^ 
( Kkms"-'^)"'^, or 3 - 7 times higher than in the Milky Way. 
We may compare this result to that obtained by Cohen et 
al. (1988). Comparing CO luminosities L(CO) to velocity 
width Au, they conclude that on average Xlmc = 6Xg. 
They adopt Xq = 0.28 x lO^i cm^^ (Kkms-^)"! result- 
ing in a value for Xlmc 30% higher than ours. However, 
if Xg = 0.20 X 10^1 cm-2 (Kkms"!)"! (Bloemen 1989), 



their result is identical to ours. Satisfactory as this may 
seem, the situation is more complex. 

First, Cohen et al. (1988) estimated their value of X 
from the mean ratio L{CO)g/ L{CO)i,mc, but their Fig. 
2 shows this ratio to be a fimction of Av. At the smallest 
velocity widths, their implied Xlmc is only 3 Xq, but at 
the largest widths it is IOXq. The figure exhibits a large 
scatter around the mean, covering an equivalent range in 
XhMC of 2 - 20 Xq. Part of this scatter is undoubtedly due 
to the low signal-to- noise ratio of the CO measurements, 
but part of it is real. For instance, Garay et al. (1993) 
studied seven CO clouds in the 30 Doradus halo and found 
those to have L{CO)g / L{CO)uot ratios implying a much 
higher value ^Dor = 20 Xq than the mean found by Cohen 
et al. (1988). We note that the values tabulated in our 
Table 2 also define a large range in X, from close to the 
Galactic value to more than an order of magnitude higher. 
For 30 Doradus itself we derive an even higher X value. 
We will discuss this variation in X in sect. 3.2. 



Table 4. Comparison of X Values 



Galaxy 


10^^ cm" 
This Paper 


X 

VKkms-i 
Literature 


References 


LMC 


1.3 


1.7; 3.9 


1; 2 


SMC 


12 


6 


3 


NGC 55 


3 


6 


4 


NGC 1569 


16 


5 


5 


NGC 4214 


3 


2; 1 


6; 7 


NGC 4449 


0.8 


1 


8 


NGC 6822 


5 


0.6-3 


9 


NGC 604 


2.2 


3; 0.35; 3.5 


10; 11, 12; 13 


NGC 595 


1.2 


1.4; 0.6 


10; 12 


NGC 5461 


1.1 


2 


14 



References for other X values: 1. Cohen et al. (1988); 2. Garay 
et al. (1993); 3. Rubio et al. (1991); 4. Dettmar & Heithausen 
(1989); 5. Grcvc ct al. (1996); 6. Thronson ct al. (1988); 7. 
Becker et al. (1995); 8. Estimated from Ohta et al. 1993, af- 
ter correction for main-beam efficiency; 9. Wilson (1994); 10. 
Estimated from Blitz (1985); 11. Viallefond et al. (1992); 12. 
Wilson & Scoville (1992); 13. Israel et al. (1990); 14. Estimated 
from Blitz et al. (1981). 

Second, we differ in individual cases, even though the 
mean values agree. For instance, Cohen et al. (1988) ob- 
tain very high H2 and virial masses (and their Fig. 2 im- 
plies a high value for X) in cloud 35. In this cloud, we find a 
low Xvalue. The result by Cohen et al. (1988) follows from 
the high Av = 28 km s~^ they found for this cloud. They 
list similarly large velocity widths for e.g. clouds 13, 19 and 
23. Yet, the higher resolution SEST survey yields velocity 
widths typically a factor of two less (Israel et al. 1993; 
Kutner et al., 1997). At least in the case of cloud 35, the 
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anomalously high velocity width appears to be the result 

of CO clouds at two distinct velocities blended together in 
the large beam used by Cohen et al. (1988). Reduction of 
the large velocity widths to the more modest SEST val- 
ues, yields X values in much better agreement with those 
in Table 2. Similar comments apply to the SMC results by 
Rubio et al. (1991), except that here we find a Zar^rer value 
of X, although both results have significant uncertainties 
associated with them. It is nevertheless clear that Xsmc 
is much higher than either Xlmc or Xq. 

For the objects listed in Table 3, X estimates can be 
obtained from the literature (Table 4). These are mostly 
rough estimates based on comparison of CO luminosities 
and virial masses, and are very uncertain. Nevertheless, 
we see reasonably good agreement in Table 4 for NGC 55, 
NGC 4214, NGC 4449 and for NGC 5461. There is also 
good agreement for NGC 595 and NGC 604 if we disregard 
the estimates derived from the high resolution observa- 
tions which apply to individual cloud components rather 
than whole complexes. It has been noted before (Rubio 
et al. 1993; Verter & Hodge, 1995) that such observations 
always yield X values much lower than the global values 
derived from observations covering the whole complex (see 
als sect. 1.2). Our results for NGC 6822, and especially for 
NGC 1569, are higher than the other published estimates. 

3.2. Dependence of X on environment 
3.2.1. Dependence on radiation field 

Compared to the Milky Way, the galaxies studied here 

have lower metallicities, and are found to have higher X 
ratios. This is not a new result: several authors have sug- 
gested a more or less inverse linear dependence between X 
and metallicity as measured by the oxygen abundance (c/. 
Dettmar & Heithausen, 1989; Rubio et al. 1991; Verter & 
Hodge 1995; Arimoto 1996, Sakamoto 1996). Arnault et 
al. (1988) found the ratio of CO to HI emission in a sample 
of 19 late-type galaxies to vary as roughly as [0]/[H]^-^. 
Sage et al (1992) failed to reproduce such a relation, but 
our sample suggests the CO to HI ratio to be propor- 
tional to [0]/[H]^-^, within the errors identical to the re- 
sult obtained by Arnault et al. (1988). The latter conclude 
that CO must be deficient relative to H2 as a function of 
metallicity, but could not determine a functional depen- 
dence for the CO / H2 ratio, hence X, on metallicity. Their 
result nevertheless suggests a dependence with a coeffi- 
cient larger than unity, as there is no reason to assume 
vastly different H2/HI ratios in low- metallicity galaxies. 
The present data provide an excellent basis to pursue this 
question, as they have been obtained in a consistent man- 
ner. Much of the previously published discussions were 
based on a compilation of data (notably X-valuc;s) from 
different sources, and obtained in different manners. 

An underabundance of CO with respect to hydrogen 
is expected to result from low carbon and oxygen abun- 



Fig. 3. Dependence of X on radiation field as represented by 
the ratio ctfir/A'h. Left: LMC; right SMC. Linear regression 
lines are plotted. 



dances. It will be enhanced by photodissociation of CO 
since low metal abundances imply both lower selfshield- 
ing and lower dust shielding of CO against the ambient 
radiation field. The LMC and SMC samples allow us to 
first investigate the effect of the radiation field. In Fig. 3 
we have plotted the ratio N(H2)/I(C0) = X as a func- 
tion of the energy available per H nucleon, represented by 
the quantity ctfir/A'h- Straight lines indicate the linear 
regression lines. In case of the LMC, we did not include 
30 Doradus in determining the regression line, because its 
very high values would dominate the result. Nevertheless, 
extrapolation of the regression line to the ctfir/ iVn value 
of 30 Doradus predicts it to have X = 7 x 10^°cm~^ cm~^ 
(Kkms~^)~^, or 85% of the value derived directly in Ta- 
ble 1. The SMC sample, although much smaller, shows a 
similar behaviour. Further analysis suggests that the de- 
pendence of X on radiation field is indeed very close to 
linear: X oc ( o-pm/ A/h)°-^=^°-^. 

The increase of X, i.e. the decrease of CO relative to 
N(H2), as more energy per nucleon is available is the re- 
sult of two processes, as is illustrated by the SEST results 
obtained for the LMC. High- resolution maps (linear beam- 
size corresponding to 10 pc) were obtained of clouds 35/36 
(south of 30 Doradus - Kutner et al. 1997), cloud 6 (Nil - 
see Israel & dc Graauw 1991) and cloud 32 (30 Doradus - 
Johansson et al., in preparation). The map of cloud 35/36 
shows numerous clumps embedded in extended interclump 
gas; the average peak-to-diffuse CO contrast ratio is about 
3. Bright clumps (i.e. those having a CO strength of more 
than 5 Kkms~^ per SEST beam area) are numerous, but 
contribute only about a quarter to a third of the CO lu- 
minosity of the whole complex. This is similar to Galac- 
tic giant molecular clouds, where e.g. Heyer et al. (1996) 
find most of the molecular mass to be in extended, low 



F.P. Israel: H2 and its relation to CO in the LMC and other magellanic irregular galaxies 



9 



column-density regions. Cloud 6 is embedded in a four 

times stronger radiation field (Tabic 1) and contains a 
very similar number of clumps per unit area. Two thirds 
of these are bright with I (CO) > 5 Kkms-^ per SEST 
beam, but cloud 6 lacks the extended diffuse gas seen in 
cloud 35/36. In this complex, the contrast ratio is of order 
ten. Apparently, the fourfold increase in radiation density 
has resulted in the removal of virtually all the low column 
density CO gas. As the high column density hydrogen gas 
will be practically unaffected, the value of X has increased 
in cloud 6 more or less commensurate with the increase in 
radiation density and CO removal. 

Cloud 32 experiences a radiation density a factor of 
three over that in cloud 6, i.e. over an order of magnitude 
more than that of Cloud 35/36. There is no trace of in- 
terclump gas. The number of clumps per unit area is still 
very similar, but the fraction of bright clumps is only a 
third, down by more a factor of two from Cloud 6. Many 
of the weaker clumps are hardly discernible. We conclude 
that the further increase in radiation density in cloud 32 
is eroding even the dense CO clumps that are surviving 
reasonably well in cloud 6, resulting in a further decrease 
of the I(CO)/ N(H2) ratio, i.e. a further increase in X. 

It is of interest that the low-cxcitation clouds south 
of N 159 have X values only a few times lO^'^cm"^ cm~^ 
rather similar to the canonical value of X in the Milky 
Way. In these clouds, the lack of dissociating radiation 
apparently allows the CO to fill most of the H2 volume, 
nothwithstanding the lower CO abundance. 

3.2.2. Dependence on metallicity 

The ratio of X to crFiR/.ArH for the objects in Tables 1, 
2 and 3 can now be compared to the mctallicities given 
in Table 5. We have included the Milky Way by using 
the Y and a-pm/ -^h values tabulated as R2 tot R5 by 
Bloemcn et al. (1990 - their Table 3) and the mctallici- 
ties from Shaver et al. (1983). We have determined linear 
regressions with and without the Milky Way points, and 
for X-a-pin/ Nyi dependences with exponents 0.9 and 1.0. 
We find: 

log X= 0.9±0.1 log(-^) - 3.5±0.2 logjgj + 34.6±2.2 
(n = 10; regression coefficient 1^ = 0.78) (3a) 

Inclusion of the Milky Way points changes this to: 

log X= 0.9±0.1 log(^) - 3.2±0.1 logjgj + 34.3±3.1 
(n = 14; = 0.85) (3b) 

This is illustrated in Fig. 4a, which shows the depen- 
dence on [0]/[H] for a linear (exponent 1.0) dependence 
oi X on (Tfir/ ^Vh. As is clear from Fig. 4a, there is almost 
perfect agreement of the Milky Way points with the re- 
lation determined for the sample galaxies alone. We have 
also empirically determined the dependence of X on metal- 
licity, ignoring any dependence on radiation density. This 
yields (see also Fig. 4): 



Fig. 4. Dependence of X on metallicity. a. At left: the ratio 
of X to ctfir/ A'h as a function metallicity [0]/[H]. The linear 
regression is drawn as a solid line. Only one line is drawn, as 
inclusion of the Milky Way points does not perceptibly change 
the result, b. At right: X as a function of metallicity [0]/[H] 
regardless of ambient radiation field. Global Milky Way points 
arc indicated corresponding to ' V from Bloemen at al.(1990 - 
filled circle) and the more commonly used 'JT from Bloemen 
et al. (1986 - open circle). Regression lines are marked for 
the sample galaxies only, and for the sample galaxies with the 
addition of the Galactic ' F point (steeper line) . 

log X = -2.7±0.3 logjgj + 11.6±1.0 

(n = 10; = 0.90) (4) 

Here, the nominal global Milky Way y-point (Bloemen 
et al. 1990) is low compared to the relation defined by the 
sample galaxies, while the commonly used value X = 2.3 
X 10^° cm~^( Kkms~^)^^ provides a very good fit. 

The dependence of X on [0]/[H] alone, ignoring 
o'fir/ .^h effects, found here is significantly steeper than 
the result log X oc 1.5 log [0]/[H] foimd by Sakamoto 
(1996) from modelling radiative transfer and excitation of 
CO in clumpy molecular clouds. However, that result did 
not take into account the full effects of photodissociation 
of CO especially on the interclump gas. 

With respect to steep dependences on metallicity, we 
note that Garnctt et al. (1995) have shown that [C]/[0] 
oc [O]/[H]0-43±0 09, so that [C]/[H] should be proportional 
roughly to [0]/[H]^-^. If the CO abundance is solely deter- 
mined by the carbon abundance, [CO]/[H] likewise will be 
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Table 5. Adopted Abundances 



Galaxy 


[0]/[Hl 

L J / L J 

10"" 


References 






LMC 


2.60±0.40 


1, 2, 3, 4 


SMC 


1.05±0.20 


1, 2, 3, 4 


NGC 55 


2.15±0.25 


5, 6, 7 


NGC 156<) 


i Qfi+n 1 5 


2 5 


NGC 4214 


2.05±0.15 


2, 7, 8 


NGC 4449 


3.0±1.0 


2, 9, 10 


NGC 6822 


1.60±0.15 


2, 5, 11, 12, 13 


NGC 604 


2.40±0.30 


14, 15, 16 


NGC 595 


2.40±0.30 


14, 16 


NGC 5461 


3.30±0.40 


9, 17, 18 



References for abundances: 1. Dufour (1984) and references 
cited; 2. Skillman et al. (1989a) and references cited; 3. Camp- 
bell (1992); 4. Pagel et al. (1992); 5. Talent (1980); 6. Web- 
ster & Smith 1983; 7. Stasinska ct al. (1986); 8. Kobuhiicky & 
Skillman 1996; 9. McCall (1982); 10. Hunter et al. (1982) 11. 
Skillman et al (1989b); 12. Kinman et al. (1979) 13. Pagel et 
al. (1980); 14. Kwitter & AUer (1981); 15. Diaz et al. (1987); 16 
Vilchez et al. (1988); 17. Rayo et al. (1982); 18. Evans (1986). 



proportional to [0]/[H]^-^; if the oxygen abundance plays a 
role this may increase to [0]/[H]^-^. The strength of the ra- 
diation field experienced by CO is proportional to the pho- 
ton flux diluted by dust extinction. To first order, we may 
equate the decrease in dust shielding with the decrease 
in dust abundance. As the dust-to-gas ratio in galaxies is 
about proportional to [0]/[H] (see e.g. Issa et al. 1990), 
we expect photodissociation alone to gain in importance 
as roughly [0]/[H]~^ when metallicity decreases. Because 
the effects of photo-dissociation are highly non-linear, and 
depend critically on the balance between ambient radia- 
tion field and local column densities, a more quantitative 
estimate of the effect of metallicity can only be obtained by 
detailed modelling, which should treat photodissociation 
more rigorously and take into account the structure and 
column density distribution of the molecular clouds ex- 
periencing the dissociating radiation (c/. Maloney 1990a). 
This is beyond the scope of this paper. 

3.3. H2 masses 

The results given in Tables 1 and 2 show the presence 
of significant amounts of molecular hydrogen in both the 
LMC and the SMC. At the (CO-selected) positions sam- 
pled, H2 locally dominates the interstellar gas. Table 3 
suggests that this is also true in the other galaxies. 

Can we extrapolate the results obtained so far to esti- 
mate the total amount of H2 in the sample galaxies? Ap- 
plication of the mean X value from Table 1 to the LMC 
CO results obtained by Cohen et al. (1988) yield M{H2) 



= 0.8 X 10^ Mq. A more detailed treatment, applying the 

individual N(H2 ) values from Table 1 to the cloud complex 
sizes given by Cohen et al. (1988) and correcting for their 
CO sources not included in our Table 1, yields a higher 
value M{H2) = 1.2 x 10* M©. This result is very simi- 
lar to that obtained by Cohen et al. (1988) (but see sect. 
3.1). It corresponds to a global mass ratio of molecular- to- 
atomic hydrogen of 0.2, much lower than the mass ratio 
of 1.9 found for the individual CO clouds in Table 1. It 
implies that in the LMC, about 12% of all the interstellar 
gas, including helium, is in molecular form. 

We have also applied eqn. (2) to the total far-infrared 
emission of the LMC, yielding a total hydrogen mass 
= 1.1 X 10* Mq, much less than the observed total HI 
mass M(HI) = 5.4 x 10* Mq (McGee & Milton 1966). 
Apparently, the LMC contains a significant fraction of 
relatively cool dust not significantly contributing to the 
total infrared luminosity. However, we may still estimate 
the total amount of H2 associated with the warm dust by 
assuming the mean H2 /HI mass ratio from Table 1 to ap- 
ply to all sources of warm H2: we then find M(i?2) = 0.7 
X 10* Mq. This rougher method thus underestimates the 
actual amount of H2 by about 30%. 

Following the same procedures for the SMC, we find 
M(H2) = 0.50±0.05 X 10* Mq, implying M{H2)/M(HI) 
= 0.1, or 7% of all interstellar gas in molecular form. This 
is a lower limit, since the CO observations by Rubio et al. 
(1991) cover only part of the SMC. Indeed, extrapolation 
to the total infrared luminosity of the SMC yields M^ = 
1.75 X 10* Mq which nevertheless still falls short of the 
total HI mass M(HI) = 5 x 10* Mq (Hindman 1967). 
Applying the mean H2/HI mass ratio from Table 2, we 
obtain M(i?2) = 1-0 x 10* Mq. 



Table 6. Global H2 mass estimates 



Galaxy 


M(HI) 


10* Mq 


A/(H2) 


M(H2) 
M(HI) 


M(H2) 
A^gas 


LMC 


5.4 


1.1 


1.0 


0.19 


0.12 


SMC 


5.0 


1.8 


0.75 


0.20 


0.12 


NGC 55 


18.6 


4.8 


2.9 


0.16 


0.10 


NGC 1569 


1.4 


0.65 


0.5 


0.35 


0.20 


NGC 4214 


11.2 


1.4 


1.0 


0.09 


0.06 


NGC 4449 


55 


21 


8.7 


0.16 


0.10 


NGC 6822 


1.5 


0.5 


0.4 


0.27 


0.16 


Mean 
Milky Way 


48 




12 


0.20 
±0.04 
0.25 


0.12 
±0.02 
0.15 



The galaxies listed in Table 3 were not sampled exten- 
sively in CO, so that we can only extrapolate from the to- 
tal infrared luminosity. However, the example of the LMC 



F.P. Israel: H2 and its relation to CO in the LMC and other magellanic irregular galaxies 



11 



suggests that this extrapolation is accurate to about 30%. 

The results arc given in Table 6, which also includes the 
global results for the Milky Way given by Bloemen et al. 
(1990). Table 6 shows that molecular hydrogen, although 
dominating the interstellar medium near star formation 
regions, occurs much less predominantly in irregular dwarf 
galaxies as a whole. Globally, the total mass of atomic hy- 
drogen is typically five times higher than that of molecular 
hydrogen. On average, about one eigth of the interstellar 
gas mass is in molecular form. These fractions are surpris- 
ingly close to those of the Milky Way Galaxy as a whole, 
where the fraction of molecular gas reaches a peak of 0.46 
in the 'molecular ring' at R = 4 - 8 kpc, and drops to 
0.11 in the outer galaxy (Bloemen et al. 1990). If Xcai is 
somewhat lower than assumed by Bloemen et al., as has 
been suggested by e.g. Bhat et al. (1986), the similarity 
of the Milky Way and irregular dwarf mean ratio is even 
more striking. Our result does not reproduce the apparent 
dependence of global molecular gas fraction on metallicity 
discussed by Tosi & Di'az (1990) and Vila-Costas & Ed- 
munds (1992). We note that the latter expressed doubts 
on the physical significance of that result, and suggested 
that it might be an artifact of the CO to H2 conversion 
used. Our result implies that this is indeed the case. 

The molecular hydrogen fraction in this admittedly 
small sample appears to be imcorrclatcd with metallic- 
ity, hence presumably dust-to-gas ratio. This is somewhat 
surprising as H2 molecules form on dust grain surfaces, so 
that one would expect less H2 in low metallicity environ- 
ments poor in dust grains. Our result may thus indicate 
that the formation of H2 is so efficient, that it is to first 
order independent of the dust abundance. Alternatively, 
it may reflect a selection effect. Higher metallicity envi- 
ronments provide more shielding and therefore may have 
a larger fraction of cold dust/molecular hydrogen than 
lower metallicity environments. In the presence of warm 
dust, IRAS fluxes poorly sample cold dust, so that we may 
increasingly have underestimated the total amount of H2 
for the higher metallicity galaxies. 

4. Conclusions 

1. IRAS far- infrared surface brightnesses and HI column 

densities arc used to indepcntly estimate H2 column 
densities towards CO clouds observed in the LMC and 
SMC. Generally, in these clouds H2 mass surface den- 
sities exceed those of HI by a factor of about 1.5 on 
average. This is in contrast to the global H2 to HI mass 
ratios which are of the order of 20±10 %. 

2. By combining the newly derived H2 column densities 
with published CO intensities, the CO to H2 conver- 
sion factors X are determined to bo Xlmc = 1.3±0.2 
X 10^^ molecules cm~^ (Kkms^^)^^ and Xsmc = 
12±2 X 10^1 molecules cm'^ (Kkms-i)-^ The global 
mass of (warm) molecular hydrogen is estimated to be 
M(H2) = l.OibO.3 X 10^ Mq for both LMC and SMC. 



3. On average somewhat higher molecular to atomic hy- 
drogen mass surface densities are found in the irregular 
dwarf galaxies NGC 55, NGC 1569, NGC 4214, NGC 
4449 and NGC 6822, as well as in the extragalactic 
HII region complexes NGC 604, NGC 595, both in M 
33, and NGC 5461 in M 101. The X-values derived for 
the HII regions and NGC 4449 are comparable to that 
of the LMC, while the X-values derived for NGC 55, 
NGC 4214 and NGC 6822 are typically two to four 
times higher; NGC 1569 has a very high value compa- 
rable to that of the SMC. 

4. Analysis suggests that the CO to H2 conversion factor 
X is linearly dependent on the strength of the ambi- 
ent radiation field per nucleon, and inversely depen- 
dent on a steep function of metallicity [0]/[H]: log X 
= 0.9±0.1 log-^ - 3.5±0.2 log|^ + 34.6±2.2. If the 
dependence on radiation field is neglected, the relation 
log X = -2.7±0.3 log|g| + 11.6±1.0 also fits the data. 
Similarly derived Milky Way values fit these same re- 
lations. They arc interpreted as the result of selective 
photodissociation of CO under conditions of high radi- 
ation field energy densities and poor shielding and self- 
shielding in low-metallicity environments. Thus, over 
the parameter range studied, the CO content of galax- 
ies varies strongly as a function of conditions. 

5. Estimates of the global (warm) H2 to HI mass ratios 
and the (warm) H2 gas fractions yield very similar 
results for all galaxies. On average, Af(H2) = 0.20 
M(HI), and M{ H2) = 0.12 Mgas. These ratios are very 
close to the global Milky Way ratios: the global warm 
H2 fraction in irregular dwarf galaxies appears to be 
very similar to that of our Galaxy, notwithstanding the 
large differences in total mass, luminosity, metallicity 
and observed CO luminosity. 
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